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ABSTRACT 

Using Milky Way data of the new Effelsberg-Bonn H i Survey (EBHIS) and the Galactic All-Sky Survey (GASS), we present a 
revised picture of the high-velocity cloud (HVC) complex Galactic center negative (GCN). Owing to the higher angular resolution 
of these surveys compared to previous studies (e.g., the Leiden Dwingeloo Survey), we resolve complex GCN into lots of individual 
tiny clumps, that mostly have relatively broad line widths of more than 15kms"'. We do not detect a diffuse extended counterpart, 
which is unusual for an HVC complex. In total 243 clumps were identified and parameterized which allows us to statistically analyze 
the data. Cold-line components (i.e., Avf„hra < 7.5 km s"') are found in about 5% only of the identified cloudlets. Our analysis reveals 
that complex GCN is likely built up of several subpopulations that do not share a common origin. Furthermore, complex GCN might 
be a prime example for warm-gas accretion onto the Milky Way, where neutral H i clouds are not stable against interaction with the 
Milky Way gas halo and become ionized prior to accretion. 
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c^ 1. Introduction 



In recent years it has become more and more clear that the gas 
in the Milky Way (MW) halo plays an important role for under- 
standing the evolution of our host galaxy. Matter is constantly 
expelled from the Galactic disk via winds and supernovae-driven 
outflows and at the same time material is accreted to the disk. 
Furthermore, gas produced by tidal interaction or remnants from 
the r nerging history of the MW may still reside in the halo 
(e.g.. ISancisi et al.]l2008l and references therein). Today, it is 
known that the so-called intermediate- and high-velocity clouds 
and complexes represent various stages of these processes (e.g., 
iKalberla & Kerpll2009l and references therein). A prime exam- 
ple for tidally disrupted material is the Magellanic Stream (MS) 
that extends over hundreds of degrees over the sky. In many other 
cases the origin of the gas is still debated. 

One of these cases is the high-velocity cloud (HVC) complex 
Galactic center negative (GCN) which i s one of the smaller cloud 
comp lexes (in total flux and spatial size;'Wakker & v an WoerdenI 
Il99lh . Owing to the lack of suitable background sources (stars or 
QSOs), relatively little is known about its chemical composition, 
gas properties, and distance. 

There are several early studies of complex GCN. The first 
H I detecti on of an individual clo ud belonging to GCN was de- 
scribed in lSaraber & Shana (11974) . This object is located close 
to the Galactic center ((l,b) = (8°, -4°)). Several origin sce- 
narios were discussed in Saraber & Shane's work, (1) an extra- 
galactic origin (with the cloud being a galaxy), (2) that it is an 
object penetrating the MW disk, and (3) that it is ejected from 
the Galactic center. In any case, the cloud would be located at a 



distance of several kpc (using some simplifications). The same 
cloud was independently studied by Mirabel & Franco (197(5^. 
They arrive at the conclusion that the object was ejected from 
the Galactic center and exclude other potential scenarios, e.g., 
that it is an ex tension of the M agellanic Stream, which is clearly 
contradicting iGiovanellil (1 198 ID . who favored the MS scenario. 
Another GCN cloud was observed by Cohen & Mi rabel (197^ . 
Again, several origins were discussed, with an additional sce- 
nario coming into play, i.e., that the cloud is a local group object 
(but not a galaxy). 

Even the most recent analyses of the neutral hydrogen 
in GCN are based on observa t ions d ating back to about 
three decades ago (iBaiaia et al.l 119851: iHuIsbosch & Wakken 
1988) which were performed with the Instituto Argentino de 
Radioastronomia 30-m and the Dwingeloo radio telescope. Both 
surveys were not fully sampled and hence are missing informa- 
tion. Furthermore, the small dish sizes led to a poor angular res- 
olution of about 30'. In these data sets GCN appears to be a pop- 
ulation of widely scattered, small clouds at high negative veloc- 
ities with sizes between 1 and 15 deg- (Wakker & van Woerdea 
|1991|) . Most authors place GCN at a distance of more than about 
10 kpc such that the size of the complex is a t least on the same or- 
der (e.g., Saraber & Shane 1974; Wakker & van Woerdeiill 19971: 
Weineret al. 2001; Wakker 2004) though the uncertainties are 
high because no direct distance estimates are available. 

Using the X-ray shadowing technique, iKerp et all (1 19991) 
could show that in the line of sight toward the galaxy Mrk 509 
(which is located in the GCN footpr int) there is excess so ft X- 
ray emission. For the same sight line lSembach et all ( Il995h per- 
formed HST absorption-line measurements and detected highly 
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Fig. 1. Column density map of HVC complex GCN. Spurious emission is visible along the MW plane that is caused by the increased 
noise level in this direction. The column densities were computed over the velocity interval visr = -360 ... - 160 km s"' using a 
mask to suppress image artifacts caused by noise. The mask was created using a hanning filtered data cube (hanning-9) and applying 
a 2.5cr threshold. The graticule shows Galactic coordinates. 



ionized high-velocity gas. The X-ray data sug gest that there is 
collisionally ionized gas, while ISembach et aP (1995) may have 
detected the cooler portion of the collisionally ionized plasma. 

A few additional sight lines revealed highly ionized gas 
(e.g., Civ and Siiv, Wakker 2001, and references therein). 
The presence of these metals provides strong evidence that 
GCN is not of primordial origin. Some authors conclude 
that GCN represents an inflow of gas to the Galactic center 
(lMirabellll982HMi rabel & Morraslll98 4 as well as the complex 
Galactic c enter positive (GCP), containing the famous Smith 
Cloud (see iLockman et al]|2008l for a re cent study). Based o n 
Lei den/ Argentine/Bonn survey data (LAB: lKalberla et al.l2005l) . 
Ijini (2010) used several clouds in the direction of complex GCN 
and calculated a hypothetical orbit that these objects might fol- 
low concluding that GCN likely was tidally stripped of a dwarf 
galaxy. 

results obtained using 
All-Sky Survey (GASS; 



In this work we present new 
the recently complete d Galactic 



iMcClure-Griflidis et al.1 l2009t iKalberla et all I20T0I) obtained 
with the Parkes telescope and data from the ongoing EfFelsberg- 
Bonn Hi Survey (EBHIS; Winkel et al. 2010a; Kerp et al. 201 1) 
measured with the 100-m telescope at EfFelsberg. The better an- 
gular resolution compared to previous large-area H i surveys and 



the full sampling of the data allowed us to obtain a revised pic- 
ture of complex GCN. 

The paper is organized as follows. In Section|2]we describe 
the data from EBHIS and GASS that were used for the analysis. 
We compiled a catalog of clouds (Section|3]l to perform a statis- 
tical analysis, the results of which are presented in Section|4] In 
Section|5] we discuss the results and conclude with a summary 
and outlook (Section|6]l. 

2. Data 

For our analysis we used data from the EBHIS and GASS sur- 
veys. In Fig. [T] we present a column density map covering the 
region l?*" < a < 23'', -50° <6 <25°. This is the area in which 
we found Hi gas likely belonging to complex GCN. Above 
6 > -5° Fig. [T] shows EBHIS data, the remaining part is GASS 
data. We also checked the overlap regime between both surveys 
(-5° < 6 < 1°) for consistency of the intensity calibration. 
Unfortunately, above 6 > 25° EBHIS data were not yet available . 
GASS data were reduced as described in IKalberla et al.l (|2010|) . 
The resulting angular resolution is 15!6, yielding an RMS of 
57 mK per spectral resolution element (Av - 0.8 km s"', effec- 
tive velocity resolution: ^v = 1 km s~'). For EBHI S the data 
reduction scheme presented in IWinkel et al.l (1201 Oal) was used. 
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sity map of the CND region (in Galactic coordinates, bottom 
panel) and a p-v diagram (top panel). The pronounced feature is 
marked with (a). We will come back to this in Section l4~4l 
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Fig. 2. Column density map of the negative velocity part of the 
central nuclear disk (CND) of the Milky Way (integrated over 
visr = -290 ... - 160 km s"') and an integrated position-velocity 
map (-2° < b < +2°). The structure labeled with (a) marks 
a large cloud association/filament located next to the CND (in 
position-velocity space). 



EBHIS has a higher nominal noise level of < 90 mK and lower 
spectral resolution (channel separation: Av = 1.3 km s', effec- 
tive velocity resolution: 6v - 2.1 kms"') but owing to the bet- 
ter angular resolution of 10'5, the resulting column density de- 
tection limit (after angular smoothing to the Parkes beam) is 
only about \^ highefl A^J^f = 4.1 ■ 10"^ cm^^ (GASS) and 



Arlimit 
"hi 



5.9 ■ 10'**cm-2 (EBHIS, unsmoothed: 8.9 ■ 10"* cm^^) 



calculated for a Gaussian-like emission line of width 20kms" 
(FWHM) and a detection threshold of 5<j. However, one has to 
keep in mind that (spatially) unresolved sources would be af- 
fected by beam-smearing degrading the detection probability, 
which in this respect favors the EBHIS survey because it has a 
better angular resolution. Bot h data sets were correc ted for stray 
radiation using the method of iKalberla et al.l (Il980h . which was 
previously applied to the LAB survey. 

Figure[T]shows the column density map of the data integrated 
over the velocity interval visr - -360 ... - 160km s"'. Note that 
for b ~ 0° there is spurious emission visible because the spec- 
tra covering the Milky Way disk have a higher noise level and 
in some cases strong continuum sources may produce baseline 
artifacts. In the lower right of the column density map, at about 
{l,b) = (-5°,0°), there is a bright feature that can be attributed 
to the central nuclear disk (CND) of the M ilky Way. The CND 
is also prominently visible in CO data; see iDame et al.l (1200 ih . 
A closer inspection of the data cube reveals an additional feature 
that cannot directly be attributed to the CND but may somehow 
be related to the GCN clouds. In Fig.|2]we show a column den- 



' EBHIS will have similar sensitivity as GASS after the completion 
of the second coverage. 



3. Compilation of a cloud catalog 

In order to build a catalog of clouds we searched for features in 
excess of 5cr,-n,s in the original data cubes and hanning-filtered 
versions of the spectra. This spectral smoothing was performed 
with the miriad task banning using kernel widths of 3 and 9 
spectral channels to optimize the search for cold and warm line 
components, respectively. However, we detected no additional 
features in the original and hanning-3 cube that were not ap- 
parent in the hanning-9 data cube, because all cold-line compo- 
nents that were found lie on top of a detection of a warm neutral 
cloud. For each identified clump we extracted the peak spectrum 
and used the GILDAS package class to fit Gaussians to the 
line profiles. If possible, two separate Gaussians were used to 
account for cold and warm (or multiple) components. 

Based on our cloud sample we find that complex GCN 
extends from I ~ 0°...70°; b ~ -60° . . . 10° (and not as 
previously reported from / ~ 0° . . . 50°; b ~ -40° . . . 10°; 
IWakker & van Woerdenl[l99ll and references therein) and cov- 
ers a fairly wide velocity range between -350 < visr < 
-70kms"'. The latter limit is not very well defined; at veloc- 
ities of v'lsr > -200 km s' the CND feature and its environment 
overlap with the GCN. Therefore, to classify whether a clump is 
part of GCN, we chose a hard limit of vi^y < -70 km s ' with the 
additional constraint that GCN clumps must be well-separated 
in position-velocity space from the CND. For the subsequent 
analysis we converted radial velocities to the Galactic standard 
of rest (GSR) frame, which is better suited for objects located in 
the MW halo, especia lly given the la rge spatial extent of com- 
plex GCN. Note that iWakked (|2004 detected GCN clouds in 
a smaller LSR velocity range of -350 < visr < -170kms '. 
The newly identified clouds with Vki > -170kms"' lie mainly 
at high (negative) Galactic latitudes and have similar properties 
and more importantly are consistent with the GSR velocities of 
the previously known clumps. 

In total we were able to identify and parameterize 243 
clumps (48 EBHIS, 195 GASS). This is about twice the number 
of objects detected in the most recent survey of GCN using LAB 
data (Kalberla & Haud 2006), which can be attributed mainly to 
the superior angular resolution and full spatial sampling, because 
most of the newly found clouds are still not resolved even with 
the smaller beam sizes of the Effelsberg and Parkes telescopes. 

Many of the clouds are connected, forming larger struc- 
tures or filaments in the column density map; compare Fig.[T] 
However, the inspection of the data cubes shows that many of the 
components of these features are well-separated in velocity and 
the filaments are likely just the superposition of several more or 
less isolated clumps. It is important to keep in mind that the col- 
umn density map shows the integrated intensity over a huge ve- 
locity range. Almost all clumps have a very small angular diam- 
eter of less than two beam widths. In a few cases we also identify 
head-tail- o r horse-shoe-like morphologies (Meverdierks 1991| 
iBriins et aDl20 00). suggesting ram-pressure interaction. 

Interestingly, we find almost no extended diffuse gas in the 
21-cm line emission surrounding the compact clumps. Because 
the column densities of a substantial fraction of the clumps are 
high both in terms of absolute value and signal-to-noise ratio, it 
is very unlikely that this is only a detection-limit effect. Several 
of the clouds show signs of ongoing interaction but some of these 
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Fig. 3. In total 33 of the isolated clumps show head-tail like morphologies. The figure shows column density maps of 16 of them 
(Clump 209 and 218 are EBHIS detections, the remaining clouds shown are from GASS data). The remaining 17 cannot as clearly 
be identified as HT structures in the column density maps, because nearby clumps partly overlap in the projection. They were 
confirmed to likely be HT clouds by a close inspection of the data cube. Contours are [3, 6, 9, 12,16, 20, 30, 40, . . .] ■ (Tims with 
cr^m, = 5 ■ lOi^cm-2 (GASS) and cr,ms = 7 ■ 10'^ cm-^ (EBHIS). 



are part of filaments or larger structures, which makes an analy- 
sis complicated. For 33 of the clouds we identified head-tail-like 
morphologies. In Fig.[3]we show column density maps of the 16 
head-tail (HT) clouds that are more or less isolated. The black 
solid lines mark the vector between the head and tail. For the 
remaining clumps the column density maps alone do not allow 
for an unambiguous definition of the HT vector Therefore, we 
inspected the data cube for a better distinction of head, tail, and 
adjacent clumps, the latter of which are superposed on the HT 
cloud in the column density maps but are sufliciently separated 
in position-velocity space. The plot also clearly shows that even 



at the higher angular resolution of the EBHIS the clumps re- 
main mainly unresolved. The HT clouds have to be treated with 
care because the superposition of two or more unresolved small 
clumps could easily mimic head-tail structures. 

4. Statistical analysis of the cloud sample 

4. 7. General properties of GCN as a whole 

In Fig. ID we show the distribution of the column densities, ra- 
dial velocities, and line widths (FWHM). Column densities lie 
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Fig. 4. Histogram of column densities, mean radial velocities, 
and line widths (FWHM) of the cloud sample. 
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in the range between 3 
dianofA^Hi - 1-4- 10'^ cm"^. The line width histogram reveals 
that the bulk of the clumps has relatively broad profiles (me- 
dian value is ~ 21kms"') and only few sources with a cold- 
gas component were found with Av < 7.5 km s' in only 11 of 
243 cases (i.e., less than 5%). This is a very exceptional case 
compared to most other HVC complexes (e.g., Kalberla & Haud 
l2006h . especially considering the fact that we find many unre- 
solved HVCs and no evidence for a diffuse H i component in the 
21 -cm emission. While in principle the lack of cold cores might 
be a detection problem caused by beam-filling, many other HVC 
complexes — even at farther distances — do reveal a substan- 
tial fraction of cold cores. If the typical size of cold cores in 
GCN is not exceptionally small, which we consider unlikely, the 
lack of a cold compon ent has to be real . Note that our finding 
agrees well with Kalberla & HaudI (12006'), who performed a sta- 
tistical analysis of multicomponent structures in HVCs based on 
the LAB survey. Some of our sight lines show line widths that 
can hardly be explained by a pure Doppler broadening. For ex- 
ample Avfwhm <: 30kms"' converts to Tq > 2 ■ 10^ K which is 
too high for a WNM cloud. This again hints at the presence of 
substructure not resolved with the telescope beams of Parkes and 
Effelsberg. 

Interestingly, other recent studies of the MW halo using 
high-res olution instrument s report on numerous small clumps 
as well. IStanimirovic et al] (2008) detected nearly 200 of these 
objects at the tip of the Magellanic Stream in precursor obser- 
vations of th e Galactic Arec ibo L-band Feed Array H i survey 
(GALFA-HI; 'Pee ket al]l2011l) . However, the authors report on 
a (somewhat) higher fraction of cold-line components of 12%. 
It will surely be interesting to re-analyze other HVC complexes 
once the new survey data are available. 




Vgsr (km/s) 

Fig. 5. Radial velocities and spatial positions of clouds in our 
catalog. We a dditionally had access to the metal absorption 
line sample of iBen Bekhti et all (l2008h . who searched for Can 
and Nai gas in QSO spectra obtained with UVES/VLT; see 
Section l4.1.2l Nondetections are marked with black crosses. 
Can detections are marked with filled triangles (velocity is 
color-coded). No Na i absorption lines were detected. The dotted 
line shows the area that we used for the analysis corresponding 
to the map in Fig.[T] 



The distribution of clouds in the three-dimensional posi- 
tional parameter space (I, b, Vgsr) is extremely wide-spread. In 
Fig.|5]we show the position of all detected clumps (in Galactic 
coordinates) with their radial velocity color-coded. Clouds in the 
northeastern part tend to have generally higher absolute values 
of radial velocity, whereas the lowest absolute values are found 
in the south. However, a close inspection of the data shows that 
even within strongly confined patchefl (of few degrees size) the 
scatter in radial velocities is large — they are in many cases dis- 
tributed over more than 100 km s ^ 



4.1.1. Column density distribution 



IChurchill et al.l (l2003h defined the column density distribution 
function (CDD), i.e., the number of detected clouds per column 
density intervajj. Fig.|6]shows the CDD as derived from our sam- 
ple. As expected, toward lower values of A^hi the catalog be- 
comes increasingly incomplete. For the range A^hi ^ 10'^ cm"^ 
a power law was fitted to the CDD, yielding a slope of -2.1 ±0.2. 
We will not compare this value to other studies at this point, be- 
cause the beam size of the instrument is thought to have some 
impact on the results. Nevertheless, the obtained slope might 



- For example at (/, b) ~ (40°, -13°), (38°, -38°). 

^ To obtain proper normalization, each value was divided by the total 
number of clouds in the sample and the bin width of each column den- 
sity interval. The latter ensures independence of the power law slope on 
the chosen binning scheme (e.g., linear or logarithmic). 
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Fig. 7. Can absorption spectra in the direction of the quasars QSOB2225-404 and QSOJ2155-0922. The upper panels show 
the corresponding Hi spectra that were spectrally binned to improve the signal-to-noise ratio. Note that Milky Way emission was 
clipped. In both cases no H i detection was made. The dashed vertical lines mark the LSR velocities of the absorbers. 
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Fig. 6. Column density distribution function of the cloud sample 
using logarithmic binning. The vertical dotted lines mark the 5cr 
column density detection limit of the data (GASS: A^Jj™' ~ 5.2 ■ 

for Gaussian-like 
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ii9 orr,-2 (,^jj ^g described by 
a power law with a slope of -2. 1 ± 0.2. 



of the sample. The regime Ahi > 10 cm 



serve as reference for future analyses of HVCs based on EBHIS 
and GASS. 

Although it is tempting to interpret the slow drop-off above 
the nominal detection limits of the surveys as a physical ef- 
fect, one must be very careful. As stated above, the column 
density detection limits were calculated for the very particular 
case of 20 km s ' line width. For clouds with higher line widths 
the limit increases (proportional to V^i'fwhm)- More importantly, 
the detection thresholds are very sensitive to the angular size of 
spatially unresolved clumps. For example, for an object with a 
size that matches the EBHIS beam, the GASS detection limit 
is enlarged by about (15'6/10'5)^ « 2.2. From our sample we 
can safely assume that a high fraction of clouds is still below 
the angular resolution of both surveys. Consequently, the detec- 
tion probability is seriously affected by beam filling problems. 
Based on our sample it is not possible, unfortunately, to con- 
clude whether the drop-off is real or not. 



4.1.2. Metal absorption spectroscopy 

Important information on the physical properties of com- 
plex GCN comes from UV absorption-line measurements. We 
had access to the extensive metal absorption line sample of 
Ben Bek hti et al.l (12008), who searched for Can and Nai gas in 
QSO spectra obtained with UVES/VLT. In total 57 QSO lines of 
sight were available for the region of interest, but unfortunately, 
they almost exclusively trace the southern part of complex GCN. 
From these 57 sight lines 10 show Can absorption with radial 
velocities that agree well with the Hi emission; compare Fig. [5] 
Only one feature at {l,b) = (-2°, -48°) has a much higher ab- 
solute velocity compared to adjacent (Hi) clumps. The median 
Can column density is about A^hi ~ 5 ■ 10" cm"^, the maximal 
value was found at (/, /?, Vgsr) = (0°, -58°, -120 kms"') and has 
A^Hi ~ 3 ■ 10'^ cm"^. In Fig. [T] we show spectra for two of the 
sight lines. In both cases no H i detection was made. 

Remarkably, no Na i absorption lines were detected. Na i is a 
tracer for the colder dense cores of clouds, while Ca ii is usually 
found in the warmer envelopes of clumps. For comparison, in 
the complete QSO absorption line sample in about 40% of the 
cases Nai was detected associated to Can absorbers. A second 
noteworthy point is that in most cases no H i emission is seen di- 
rectly on top of the absorption features. This suggests that there 
is a lot of material below the detection limit of the H i surveys — 
metal absorption spectroscopy is a very sensitive tool to detect 
the low-column density portion of the gas. 

Using spectral data from the Goddard High R esolution 
Spectrograph (GHRS) on HST, ISembach et aP (Il999h detected 
absorption by intermediate and high ions such as C ii. Si iii, C iv, 
and Si IV in the velocity range between -200 and -300 km s"' 
toward the quasar Mrk509; (l,b) = (36.0°, -29.9°). We show 
in Fig. [8] (left panel) the absorption profiles of C ii, C iv, and 
Si IV toward Mrk 50 9, based on archival HST/STIS data (see also 
iRichter et alJl2 009'). 

The Mrk 509 sightline intersects complex GCN in a region 
that is close to a larger population of Hi clouds seen in 21- 
cm emission (see Fig.[T]i. No 21 -cm emission is seen toward 
Mrk509 itself, so that Sembach et al. (1999) classify this high- 
ion absorber toward Mrk 509 as "highly-ionized HVC". From 
our 21 -cm data it is now evident that the high-ion absorption 
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Fig. 8. Absorption profiles of C ii, C iv, and Si iv toward Mrk 509 
(left panel) and PKS 2155 (right panel), based on archival 
HST/STIS data. 



toward Mrk 509 traces the ionized component of complex GCN. 
From a p hoto-ionization model of the observed ion ratios toward 
Mrk 509, ISembach et alJ (Il999h derive characteristic gas densi- 
ties of nu ~ 10"^ cm"^, gas temperatures of T = 2 ... 3 ■ 10"* K, 
and thermal pressures of P/k = 1 ... 5 cm""* K for the ionized gas 
toward Mrk 509. The total column density of ionized hydrogen 
in the HVC is estimated to be logA^Hn = 18.5 . . . 19.7 [cm^^], 
depending on the metallicity of the gas. 

The second G CN sightline t hat ha s been studied in 
UV absorption by ISembach et alJ (Il999h is the one toward 
PKS 2155-304; (l,b) = (17.7°, -52.2°). It passes through an 
outer region of complex GCN that is completely devoid of 
Hi 21-cm emission (see Fig.[T]). Toward PKS 2155-304 high- 



velocity absorption is seen between -90 and -330 km s"; see 
Fig.[8](right panel). The intermediate- and high-ion absorption in 
complex GCN between -170 and -330 km s' is much weaker 
toward PKS 2155-304 compared to the Mrk 509 sightHne, sug- 
gesting that the ionized (and thus the total) gas mass is sub- 
stantially lower in this direction. Another prominent absorption 
feature that is possibly related to complex GCN is present near 
-130kms '. Here, the weak intermediate-and high-ion absorp- 
tion is accompanied by relatively strong absorption in Si ii and 
C II (see also Sembach et al. 1999. their Fig. 3). 

Some interesting conclusions can be drawn from these num- 
bers. For instance, the Mrk 509 data suggest that the column den- 
sity of ionized hydrogen in complex GCN is substantial (cov- 
ering a much wider spectral range than traced by 21-cm emis- 
sion), possibly even higher than the (typical) neutral gas col- 
umn density in the Hi clumps detected with EBHIS and GASS. 
This is supported by the detection of high-velocity O vi absorp- 



tion toward Mrk 50 9, which extends from 
(IWakker et al.ll2003h . 



-90 to -350 km s" 



4.2. Evidence for a tidal gas origin of GCN 

4.2.1 . Possible association witli tine Magellanic System 

N-Body simulations of the Magellanic System including drag 
forces (GN96-i-drag model; Diaz & Bekki 20 111) predict gas i n 
the region covered by G CN (fi gure 2 of ,Diaz"& Bekkill201 ID . 
Because Diaz & Bekkj (1201 ih only provide the distribution 
of GSR velocities and galacto-centric distances projected to 
(Magellanic) longitude in their paper, the authors kindly gave 




Vgsr (km/s) 

Fig. 9. Radial velocities and spatial positions of clouds in our 
catalog as in Fig.|5](filled circles). Diaz & Bekki kindly allowed 
us access to their G N96+drag simulation data (filled triangles; 
iDiaz & Bekkill2011 ). The solid line is a proposed orbit for some 
of the GCN clouds based on LAB data as calculated by Ijini 



us access to their data, which allows us to compare the three- 
dimensional (I, b, Vgsr) distribution of the simulated N-Body par- 
ticles with our sample; see Fig.|9] 

Apparently the spatial position of the long "filament" (in 
the N-S stripe) almost matches the simulated data. However, 
the simulated velocities are offset by about +100 km s"' when 
compared to the median value of the data — though the scat- 
ter along the filament is huge (~ 100 km s"' difference between 
minimal and maximal values). The large scatter itself is also 
not predicted, the simulated stream of gas has a scatter of about 
~25kms"'. 

Furthermore, the two large cloud regions in the northeast- 
ern and southeastern parts are not predicted by the simulations. 
On the other hand the simulations predict a lot of gas (with 
Vgsi ^ -175kms"') in the southwestern part of the area, where 
we do not detect any cloud in the investigated area (marked with 
the dotted line in Fig.|9]l- We also briefly checked the GASS 
data for a > 23*^ (beyond the southeastern edge of the ob- 
served area) and confirm the presence of gas belonging to the 
MS (vgsi- ~ -125kms"') as visible in the figure. Note that in 
the simulations the filament is an extension of the Leading Arm, 
while the gas in the southwest is associated to the Magellanic 
Stream. Consequently, the distances of both features differ. The 
former has distances between 14 and 22kpc (higher values in 
the northern part) with a mean of 18 kpc — which would also be 
consistent with the distance estimates given for GCN. The lat- 
ter objects lie at much farther distances of 30 to 60 kpc, but they 
should still be detectable with EBHIS and GASS. Note also that 
given the geometry and galacto-centric distances, the filament 
must have penetrated the H i disk of the MW (at a radius of about 
15 kpc), which of course introduces lots of physical interaction 
processes that likely are not well described by the simulations. 
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Fig. 10. Head-tail vectors as a function of Galactic coordinates. We have differently marked clouds where Vhead < i^'taii (because all 
clouds move with negative radial velocities, this means that the tail moves slower toward the observer than the head) and Vhead > I'taii- 
In some cases the velocity fields did not show a clear sign of a separation. The length of each vector indicates the separation of the 
head and tail but was scaled with a factor of ten for improved visualization. 



4.2.2. Direct orbit calculation 



Based on LAB data, 'J^ (1201 Oh recently used several clouds at 
velocities v„si - -170. . . - 30kms~' to compute a hypothet- 
ical orbit that the gas might follow. This orbit is also shown 
in Fig. |9] marked with the solid line (color again encodes Vgsi). 
Interestingly, the proposed orbit describes the spatial position of 
most of the clouds in the N-S stripe pretty well. However, be- 
cause the radial velocity scatter is huge along the N-S stripe, 
the orbit cannot predict the velocities of all features. This could 
be because the orbit calculation could not take into account the 
complex physical interaction processes that probably play an im- 
portant role for GCN (see also Section|5]l. 

The orbit predicts a heliocentric distance of the associated 
GCN clouds in the range of abou t 15 t o 35kpc (with farther 
distances in the south, see also Jinl l2010l their fi gure 4). While 
the low er limit is consistent with the results from lDiaz & Bekkil 
( 1201 lb . the upper limit is quite different. Moreover, the sign 
of the distance gradient is opposite. Jin (2010) conclude that a 
dwarf galaxy is likely to be the origin of the tidal stream, though 
no optical counterpart was detected along the orbit so far 

4.3. Head-tail structures as tracer for the motion ofciouds 

In Fig [To] we show the head-tail vectors as a function of Galactic 
coordinates. The head of the arrows mark the position of the 
heads of the clouds. Most of the HT vectors seem to be ran- 
domly distributed, but in the northeast (/ ~ 10°... 25°; b ~ 
-30° ... - 5°) we find nine HT vectors with coherent direc- 
tions that point southward (and one apparently points north- 
ward). Furthermore, in the northern part there are six HT vec- 
tors pointing westward. Usually the HT vectors are interpreted 



as motion of the clouds relative to their environment because the 
more diffuse tail is thought to be created by ram-pressure strip- 
ping. In this scenario one would expect the velocity of the tails 
(relative to the surrounding medium) to be lower than that of the 
heads. Therefore, we calculated the velocity gradient along the 
HT vector where possible — in some cases the superposition of 
clumps made the velocity fields unusable and in other cases no 
clear gradient was visible. We have differently marked the vec- 
tors in Fig [To] according to their velocity gradient. Clouds with 
Vhead < ^taii show the cxpcctcd vclocity difference (ram-pressure 
scenario). Generally, the absolute differences of the head and tail 
velocities are relatively small (at most 6 km s"'). 

Interestingly, although the apparent motion of HT clouds in 
the northeastern and northern regions are coherent, the observed 
velocity gradients are not consistent with the velocity differences 
between head and tail expected from a simple ram-pressure sce- 
nario. About half of the vectors have the opposite sign of the ve- 
locity gradient. While one could in principle argue that overlap- 
ping clumps might mimic head-tail structures, this should lead to 
a more random alignment of the HT vectors except if there is a 
preferential alignment of adjacent clumps. This could eventually 
be explained by gas fragmenting along its motion vector. Then 
the HT vectors would not mark the motion of the clouds but at 
least the alignment of the motion. Interferometric measurements 
of the suspicious clouds might help to solve this puzzle because 
they would allow for a better distinction of the head and tail. 
Possibly, the observed effect is even natural. iKalberla & HaudI 
(2006) analyzed the typical radial velocity offsets of cores with 
respect to their envelope (i.e., cold vs. warm line component) 
and find that both faster and slower cores have approximately 
the same frequency of occurrence. 
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Fig. 11. Line width of all identified clumps as a function of Vgsr. 
The solid line is the best-fit linear relation with a value of 0.01 ± 
0.01km s"'/(kms-i). 
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Fig. 12. Column densities of all identified clumps as a func- 
tion of Vgsr. To parameterize the upper envelope, the 90% per- 
centiles for certain velocity intervals were calculated (using vari- 
able bin widths such that each bin contains the same number of 
clouds). The solid line is the best-fit linear relation through these 
data points, showing a significant gradient of about -13 +4 ■ 
10"'cm-2/(kms-'). 



A A. Correlations between Av, A^™" and Vg 
4.4.1. Full sample 



In lWinkel et al.l ( 1201 Obi) we reported about a possible line-width 
- radial-velocity relation in precursor observations made with 
EBHIS covering parts of the northern region of complex GCN. 
Consequently, one of the key questions of this work was to con- 
firm or reject this relation with the help of the EBHIS and GASS 
data that cover the full GCN field. 

Here, using the complete sample of GCN cloudlets, we do 
not find a significant global Av-Vgsi gradient; see Fig.fTTl 

In Fig. [12] we show the column densities as a function of ra- 
dial velocity of the clouds. Apparently, the upper envelope de- 
creases toward lower absolute velocities. To parameterize this 
behavior, we calculated the 90% percentiles for certain veloc- 
ity intervals (using variable bin widths such that each bin con- 
tains the same number of clouds). A linear function was fit- 
ted to the resulting data points. A significant relation between 
these upper column densities N^^^ and Vgsi- is obtained of about 
-13±4- 10i''cm-2/(kms-i). 



■o 




80 70 60 50 40 30 20 10 -10 
I (deg) 

Fig. 13. For some of our analyses we defined subsets of the cloud 
catalog. Four regions where used. Stripe N-S, Stripe (a), and the 
northeastern and southeastern regions. For the motivation of the 
different regions see text. 

However, as can be seen in Figs.|5] and [TO] one can hardly 
speak of a single-cloud population, but there seem to be several 
distinct regions with respect to /, b, and Vgsr. 

4.4.2. Subpopulations 

The first very prominent population is located along the fil- 
ament (hereafter referred to as Stripe N-S) predicted by the 
GN96+drag model (Diaz & Bekki 2011). It manifests itself al- 
most like a "wall" in the three-dimensional parameter space be- 
cause the scatter in Vgsi- is considerable along this curve. Two 
other regions appear visually distinct, one in northeastern part 
(denoted as Region NE, having Vgsi < -120kms ') and one in 
the southeast (Region SE, Vgsi > -140kms"'); see Fig. [13] The 
northeastern region contains several HT clouds that apparently 
move southward. Both regions, NE and SE, show much less scat- 
ter in radial velocities than clouds in Stripe N-S. Although the 
mean GSR velocity in the NE and SE regions is different by 
about 75kms"' and there is little overlap, the two regions are 
connected in velocity via the N-S stripe. 

We introduce a fourth subpopulation, with objects located 
along a slice defined by the orientation of the IVC feature (a); 
compare Fig. [2] This slice is hereafter denoted as Stripe (a); see 
Fig.[T3] Five HT clouds lie within Stripe (a) and surprisingly 
have their HT vectors closely aligned with the slice direction. 

Again, we parameterized possible Av-Vgsr and A^^^-Vgir cor- 
relations using linear regression. Table[T]shows a summary of all 
relevant fit parameters for the different subsamples. 

Except for Stripe (a) we find no significant correlation of 
Av vs. Vgsr. Stripe (a) also shows high evidence for a column 
density evolution with Vgsr. Because Stripe (a) contains several 
HT clouds that are morphologically aligned with the CND fea- 
ture (a), we additionally had a look at the column densities and 
line widths of the HT clouds (marked with triangles in Fig.[l4li. 
In the line width-velocity panel the HT clouds are remarkably 
tight-correlated but we emphasize that only five data points are 
available. Hence, we performed a Monte-Carlo test where we 
randomly drew five data points out of the Stripe (a) sample. In 
only 3% of all cases we obtained a correlation with similar qual- 
ity of fit (in terms of reduced ;i^^ value), which makes it unlikely 
that the observed relation in the HT subsample of Stripe (a) is 
just chance. We also point out that the correlation of the line 
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Fig. 14. Line widths (upper panel) and column densities (lower 
panel) as a function of radial velocity for clumps along Stripe (a) 
as defined in Fig.[T3] Head-tail clouds are marked with a trian- 
gle. For the line widths we computed a linear fit through all data 
points (solid line) and through HT clouds only (dotted line), re- 
spectively. For the column densities we calculated the 90% per- 
centiles as in Fig. [12] to parameterize the upper envelope. The 
solid line is the best-fit linear relation of A^™!"' vs. Vg,, neglecting 



the outlier at v„ 
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and A?"f - 

±1 1 



Table 1. Results of the linear fits for the Av- 

relation for the different (sub)samples. The corresponding plots 

are shown in Fig.[TT][T2l and[T4l 
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t Neglecting an outlier, see Fig.ll4l 



width with the radial velocity is not related to a possible evo- 
lution of the line widths with spatial position (although there is 
some evolution of v^sr with position, as seen in Fig.|5j. To show 
this, we plot in Fig.llSJthe line widths as a function of position 
along Stripe (a), which do not show a correlation. Apparently, 
the line width-velocity relation is a genuine property indepen- 
dent of the spatial position of the clouds. 
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Fig. 15. Line widths as a function of spatial position (measured 
from the east) for clumps along Stripe (a) as defined in Fig.fTSl 



The question remains why all other subsamples do not show 
a Av-Vgsi correlation. One could argue that the complexity of 
the whole field makes it difficult to find the correct subsets of 
clouds that are physically connected. This is at least suggested 
from Fig. [To] where it appears that apart from few regions the 
dynamical motions are quite stochastic. 

Although no reliable Av-Vgsr relation could be found for all 
regions except Stripe (a), the maximal column densities vary sig- 
nificantly with Vgsi. Several of the regions (Stripe (a) and N-S 
as well as the complete sample) show a significant decrease for 
the upper envelope of the column densities with lower absolute 
values of the radial velocities. This might be a sign of ongoing 
interaction with the ambient halo material that dissolves and de- 
celerates the clouds. 

4.5. Statistical properties of isolated clouds — substructure 

In this subsection we discuss the statistical properties of clumps 
that are isolated, i.e., not part of a filament or overlapping with 
other clouds. This allows us to introduce another cloud pa- 
rameter, ffsub, which shall be an indicator for the amount of 
substructure for a clump. If a cloud were point-like (or had 
a Gaussian-like column density profile), it would produce a 
Gaussian-shaped column density profile in the data as the tele- 
scopes performs a convolution of the true brightness distribution 
with the telescope beam. Fitting a Gaussian to the data should 
then lead to a residual that exhibits the same noise level as the 
surrounding (emission-free) regions. On the other hand, when 
substructure is present, a Gaussian would only provide a poor fit 
to the column densities, hence the residual noise level should be 
higher than before. We define 



ffsub = 



(Ti-njs (fit residual) 
cTji^s (surrounding) 



(1) 



as the residual noise level after fitting 2D Gaussians to 
the clumps relative to the noise of emission free regions. 
Unfortunately, owing to the limited angular resolution of the 
single-dish data, one can treat the amount of substructure only 
in this statistical way. 

In Fig.[T6]we show two examples: a clump that is fairly well- 
fitted by a single Gaussian (top panels) and another clump that 
has enhanced noise in the residual (bottom panels). In the former 
case we obtain asub = 1-1, while for the latter ojsub = 2.5. Most 
of the isolated clouds in our catalog are located in the southern 
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Fig. 16. For two isolated clumps a two-dimensional Gaussian 
was fitted to the column density profile (left panels). The resid- 
ual (right panels) shows diff'erent behavior, for the upper row the 
noise level is very similar to an emission-free region in the sur- 
rounding, the lower row exhibits remaining features increasing 
the measured noise value. The residual RMS can be used as an 
estimator to the amount of substructure within a clump. Contour 
lines in the left panels start at 3o-i-ms and have a spacing of 3o-isas- 
Contour lines in the right panels are spaced with Icrrms, solid 
contours mark positive values, dashed contours are negative, the 
thick solid lines are zero flux. 
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Fig. 17. Substructure estimator Csub for all isolated clumps as a 
function of Vgsr. There appear to be two independent populations, 
which we separated using the dotted line. For the circles a linear 
function was fitted to the data, showing a gradient of 4 + 1 ■ 
10-^(kms"')-i. 



hemisphere. Therefore, we restricted the following analysis to 
the GASS data, because it is not clear in which way the different 
beam sizes and noise levels of EBHIS and GASS would intro- 
duce systematic effects. In total the subset of isolated clumps 
contains 63 objects. 

In Fig.[T7]we show asub as a function of Vgsr. There seem to 
be two different populations, which we separated for some of 
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Fig. 18. Positions of isolated clumps in Galactic coordinates. 
Triangles and circles mark the two populations defined in 
Fig-El The substructure estimate asub is color-coded. For the 
population marked with circles we find a gradient of cKjub in N- 
S direction with higher values in the southern part of complex 
GCN. 





30 


CO 

E 


25 




20 










a 


lb 


F 




S 


10 


5 






5 








- m = -0.07?t 0.03 
n = -6 ± 6 

- RMS = 2.7 


f(x)=mx+n 
m = -0.03 ± 0.01 
, n = 11 ±2 
RMS = 5.0 

• 
• 


>• ^^^r^~ 


• 




.. ' 


■'^~~~~\ 


>^ 



u 



-200 



-150 



"gsr 



-100 
(km/s) 



-50 



Fig. 19. Quantity 



Avfwhm 



as a function of Vgsr for isolated clumps. 



Triangles and circles mark the two populations defined in 
Fig-El The substructure estimate asub is color-coded. For both 
populations linear functions were fitted to the data points. A bi- 
furcation is observed. 



the following analyses. For clumps below the dotted Une (blue 
circles) a correlation of asub and Vgsr is found: 5(asub-Vgsr) = 
4 + 1 • 10"^(kms"^)"^ In Fig.[T8] asub is plotted vs. Galactic 
coordinates. Triangles and circles again mark the different popu- 
lations as defined in Fig. El There appears to be no clear spatial 
separation between both subsets. For one of the subsets (circles) 
one can see a trend in N-S direction of asub increasing toward 
the south. While one could argue that our method of defining 
the amount of substructure might be subject to systematic ef- 
fects and uncertainties, it is remarkable that this general trend is 
observed, which can hardly be explained by coincidence. 

There is no appreciable trend of asub with the line width of 
the clouds but the correlation of asub and A^h i although the scatter 
is large (not shown here). For comparison we also computed the 
line width-radial velocity relation for the subsample of isolated 
clouds, which does not exhibit any clear correlation. However, 
based on the idea that the amount of substructure might increase 
the effectively measured line width, we display in Fig. [191 the 
quantity 



™'™ as a function of Vosi- and observe, first, a bifur- 
cation of both populations (circles and triangles) and, second, 
a correlation of this "effective" line width with radial velocity, 
which is less scattered about the linear fit than before, even so 
we did not divide the catalog into different subsamples based on 
their spatial position. 
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Note that the observed relation between a^ub and A^hi can 
be explained in a natural way because two similarly distorted 
objects with different amplitudes would not necessarily lead to 
similar values of a^ub because the residual scatter is expected to 
be higher if the Gaussian has a larger amplitude. However, in 
that case similar results as shown in Fig. [19] should be obtained 
when using ^''™'" instead of ^''™"" , but this is not the case. 
Also, the two different populations are not as clearly separated 
in the A^Hi-v'gsr space as in o'sub-Vgsr- It seems that the substruc- 
ture estimator provides a more genuine property than the column 
density despite their apparent correlation. 



5. Discussion 

The observed properties of clouds in our sample show that GCN 
is quite different compared to other large HVC complexes. Not 
only the small spatial size of most of the clumps and the low frac- 
tion of cold-line components are rather untypical, but also the 
missing diffuse extended H i emission, which is detected in most 
other HVC complexes (e.g., Kalberla & Haud 2006, using the 
less sensitive LAB survey data). This is especially true because 
the absorption line data indicate the presence of a (low-column 
density) diffuse gas phase. Furthermore, we find many signs for 
collisionally or photo-ionized gas in the region. Several head- 
tail features were observed, and we find a decrease of the maxi- 
mal H I column densities toward lower absolute radial velocities. 
This makes it likely that GCN clouds are subject to interaction 
processes with an ambient medium (i.e., the Milky Way halo). 

Unfortunately, the large distance uncertainty does not allow 
us to estimate important physical parameters like temperature or 
pressure with sufficient accuracy to rule out or accept the various 
possible interaction processes directly, i.e. Rayleigh-Taylor- and 
Kelvin-Helmholtz instabilities, ram-pressure, or shocks (which 
could also lead to fragmentation). However, we will discuss be- 
low some of the mechanisms under the assumption that clouds 
in the N-S stripe are indeed r elated to the L e ading Arm, as sug- 
gested by the simulations of iDiaz & Bekkil (1201 ll) . which con- 
strains the distance to 14 to 22kpc. This is at least consistent 
with earlier (rough) distan ce estimates and partially with the or- 
bit proposed by Jinl (l2010l) . 

In this case, for a beam-sized (unresolved) cloudlet at 1 8 kpc 
distance with a typical column density of 10'^ cm"^ we estimate 
the volume density to be on the order of Mc - 0.1 cm"^. The 
size of such an object would be about 50 pc . Furthermore, us- 
ing the MW model of iKalberlaet al.l (12007 *). we can obtain an 
estimate of the pressure and density of the surrounding halo gas 
along the N-S filament: P/ks ~ 10^ Kcm""* and Wh S 10"^ cm"^ 
close to the MW disk and P/ks < 10- . . . 10^ Kcm"^ and n^ ~ 
10 "^ . . . 10"^ cm"^ at higher Galactic latitudes. Typical masses 
of the larger structures, e.g., at {l,b) - (37°, -38°), would be 
few times 10^ Mq, the smaller clouds could have a few lO"* Mq. 

In their analy s is of the clouds at the tip of the MS, 
IStanimirovic et al.1 (|2008|) estimate time and length scales for 
the physical processes important for interacting clouds. While 
at a first glance our objects have similar properties (i.e., density, 
size, and line widths), the surrounding medium makes a big dif- 
ference. Particularly, we observe signs of ongoing ram-pressure 
interaction (the head-tail structures), which were absent in the 
MS sample. This is not surprising, given the high densities of the 
halo material at smaller distances compared to the MS. Hence, 
in our case the Kelvin-Helmholtz time scales are larger (about 
2 Gyr), while thermal instabilities in principle should have simi- 
lar time scales of a few tens of Myr. 



iHeitsch & PutmaiJ (l2009l) used simulations to study the sta- 
bility of neutral clouds against ram-pressure for different scenar- 
ios (wind-tunnel and free-fall in the MW potential). While both 
cases are not exactly applicable for the GCN clouds, some of 
the authors' basic conclusions might be of interest for us. First, 
if the initial masses are on the order of 10^ to 10* M©, the time 
scale for complete disruption is on the order of 50 to at most 
100 Myr (cloud properties in our sample are similar to their sim- 
ulations and the environmental densities are likely even higher). 
Although cooling times (and hence the time scale for thermal 
instabilities) are likely shorter or comparable, the wind-tunnel 
simulations show no significant cooling of the clouds. Because 
the clouds loose material through the ram-pressure interaction, 
cooling is even less effective. This could explain the small frac- 
tion of cold objects in our sample. Note that the ablated mate- 
rial might perhaps form smaller fragments, but these would have 
much smaller masses than the main cloud. A second important 
point is that the fragments are most likely not stable (because 
of the less efficient cooling) and will eventually evaporate. We 
conclude that ram-pressure interaction can most likely explain 
the observed cloud properties, though simulations better suited 
to the GCN scenario would be appreciated. 

Under these considerations, complex GCN might be a prime 
example for warm gas accretion onto the Milky Way, where neu- 
tral H I clouds or a stream of gas are encountering the Milky Way 
gas halo. Under certain physical conditions clouds are no longer 
stable against interaction with the ambient medium, might dis- 
solve and become ionized prior to acc retion. This scenario is a lso 
suggested by recent simulations (e.g.. lBland-Hawthornl2009l) . In 
that case the GCN objects would just be in a transient state. This 
would match our observed N^'^'^-v„sr relation. The lack of Na i, 
which is a tracer for the colder cores of typical clouds, is another 
sign that the typical life times of the clouds might not be long 
enough to allow for substantial cooling of the cores — which 
explains the under-abundance of cold clouds in our sample com- 
pared to most other HVC complexes. 

If the hypothesis is true, clouds would not be visible any- 
more, which immediately raises the question of how much ion- 
ized gas mass is contained in complex GCN. Because it is hard 
to observe, it is still unclear today whether the (global) ionized 
fraction of accreted material is able to explain the observed con- 
stant star-formation rate of t he MW. Metal absorption Une stud- 
ies (using large samples; see iBen Bekhti et ani2008l: IShull et al.l 
2011) of the HVC sky will hopefully shed light on this issue. 
One should also keep in mind that Can was detected in sight 
lines without H i emission, hence, there definitely is low-column 
density material. 

We do not only observe a general relation between the up- 
per envelope of the column densities and the radial velocity, but 
this relation is also significantly found for two subregions (Stripe 
(a) and N-S). Although the other two subsamples (Region NE 
and SE) do not individually show a significant relation, all of the 
slopes are negative, which increases the overall significance. The 
correlation for Stripe (a) and N-S is surprisingly tight consider- 
ing the huge scatter in the observed radial velocities in these 
regions, while for the regions NE and SE, which have much 
less scatter in radial velocities, the relation is much less con- 
strainted. Furthermore, given the observed (non-)alignments of 
the HT vectors in each of the regions, we conclude that clouds 
in the four subsamples may have different origins. 

One striking result is the detection of five HT clouds in 
Stripe (a) — the definition of which was based only on the IVC 
feature near the CND — pointing eastward. This is compati- 
ble with the observed A^™f - 

rl 1 



relation under the hypothesis 
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that clouds with lower velocities, i.e., which were decelerated 
for a longer period of time, have already lost more gas through 
ram pressure. Here, we cannot prove a causal relationship, but it 
could be that the IVC feature is just the oldest component within 
Stripe (a), which has just encountered the MW disk. 



6. Summary and outlook 

The analysis has shown that complex GCN is indeed a very com- 
plex object consisting of more than two hundred (mostly individ- 
ual) tiny clouds, which are often unresolved even with the much 
better angular resolution of the EBHIS and GASS data compared 
to any previous study of GCN. We do not observe diffuse ex- 
tended 21 -cm emission around the isolated clouds, although the 
column density detection limit of our data would be sufficient to 
detect this gas phase at reasonable distances. Our observations 
suggest that the clouds are likely a transient phenomenon in a 
very turbulent and highly dynamic environment as expected if 
GCN is really the result of ongoing accretion. It might be that 
the clouds become ionized prior to accretion, which might ren- 
der a large fraction of the total gas mass in GCN invisible for 
the H I survey data. The fact that we detect Ca ii signatures and 
other ionized species at velocities consistent with the H i struc- 
tures supports the hypothesis that there is a significant amount 
of low-column density material around. 

The presented work is only the first step toward a deeper 
understanding of complex GCN as a whole. We are going to 
use interferometric observations for studying the morphologi- 
cally interesting cases at much better angular resolution. This 
would allow us to obtain better constraints on the true motion 
of these clouds, because currently there is often an ambiguity 
between the morphology and observed radial velocity differ- 
ence between head and tail. It would also be interesting whether 
unresolved clumps further fragment into individual cores, be- 
cause this could provide a consistency check for our substructure 
estimator. Furthermore, deep single-dish observations of some 
smaller regions could shed light on the question of the missing 
diffuse extended component. If it turns out that the latter is not 
physically present, it might be worthwhile to search for addi- 
tional absorption lines to study the ionized gas phase to estimate 
its total mass contribution. 

While this work mostly dealt with the general (statistical) 
properties of our cloud sample, we will next analyze potential 
subpopulations in more detail (Darmstadter in prep.). Based on 
our findings, we suspect that GCN is a mixture of various ob- 
jects with different origins, which one would need to separate to 
understand the complex as a whole. 
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